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ABSTRACT 

We analyze in detail the penumbral structure found in a recent radiative MHD simulation. Near 
r = 1, the simulation produces penumbral fine structure consistent with the observationally inferred 
interlocking comb structure. Fast outflows exceeding 8kms^^ are present along almost horizontal 
stretches of the magnetic field; in the outer half of the penumbra, we see opposite polarity flux 
indicating flux returning beneath the surface. The bulk of the penumbral brightness is maintained by 
small-scale motions turning over on scales shorter than the length of a typical penumbral filament. The 
resulting vertical rms velocity at t = 1 is about half of that found in the quiet Sun. Radial outflows 
in the sunspot penumbra have two components. In the uppermost few 100 km, fast outflows are 
driven primarily through the horizontal component of the Lorentz force, which is confined to narrow 
boundary layers beneath t — 1, while the contribution from horizontal pressure gradients is reduced 
in comparison to granulation as a consequence of anisotropy. The resulting Evershed flow reaches its 
peak velocity near r = 1 and falls off rapidly with height. Outflows present in deeper layers result 
primarily from a preferred ring- like alignment of convection cells surrounding the sunspot. These 
flows reach amplitudes of about 50% of the convective rms velocity rather independent of depth. A 
preference for the outflow results from a combination of Lorentz force and pressure driving. While 
the Evershed flow dominates by velocity amplitude, most of the mass flux is present in deeper layers 
and likely related to a large-scale moat flow. 

Subject headings: convection - magnetohydrodynamics - radiative transfer - sunspots - Sun: surface 
magnetism 



1. INTRODUCTION 

Since the discovery of the Evershed effect about a 
century ago (lEvershedlfTOOl . the origin of large-scale 
outflows in sunspot penumbrae has been a central ele- 
ment in observational and theoretical studies of sunspots. 
Over the past decades, advancements in ground- and 
space-based observing capabilities have revealed the 
stunning flne structure of sunspot penumbrae that is 
manifest in the intensity, magnetic fleld and velocity 
struc t ure (see, for example, recent reviews by ISolankil 
[20031: iThomas fc Ws3 [200l l2008l and high res olution 



observations by Scharmer et al. 20021: TL 



2005'; 'Rimmele & Marino"2006'; 'Ichimoto ct al. 2007a"Jl 
Langhans ct al. 2007; Scharmer ct al. 2007; Rimmela 
IFranz &: Schhchc nmaieii 120091 : IBellot Rubio et aLl 
All quantities show in the penumbra a pri- 
marily radial filamentary structure. Strong horizon- 
tal outflows take place in regions with almost hori- 
zontal fleld, embedded in a background of more ver- 
tical field - which has been referred to as "uncombed 
penumbra" ( Solan M fc Montavon 1993,') o r "inter locking- 
comb" structure (|Thomas fc Weisiil992D . The connec- 
tion between the Evershed flow and the intensity struc- 
ture is less clear. While earlier work pointed toward a 
flow preferenti ally in the dark component , more r ecently 
[Schlichcnniaicr et al.l ([TOOS) and Ichimot o et all ()2007al) 
showed that Evershed flow and intensity variations show 
a positive correlation in the inner and negative correla- 
tion in the outer penumbra. Another controversial as- 
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pect is the dep th pr oflle of the Evershed fl ow. While 
iRimmeld ()1995D and iStanchfield et al.l (|1997[ ) found the 
Evershed flow in elevated fl o w ch a nnels, more recent 
work by [Schlichcn maier et al.l ( 2004D ; IBellot Rubio et all 
([2006') and lBorrero et al.l ([2006 ) points toward a flow pri- 
marily in the deep photosphere that declines with height. 

A variety of simplified models have been pro- 
posed to e xplain the penumbral fi n e stru c ture (e.g. 
DanielsonI 119611: iMever fc Schmidt[ [19681: iGallowav " 



19751: [Thomas' '1988': 'DegenhardiJ 11989'. '199lf 'Grosse 
199ll: Wcntzcl 1992; Thomas & Montcsinos 199' 
Montcsino s fc Thomas! | 1997|: ISc hlichcn maier et al 
1998a.b, : .Spruit fc Scharmed l200"6rTScharmer fc Spruit 
Studies of idealized u iagnetoconvecti on in 



inclined magnetic field (see, e.g., IHurlburt et all [19961 
l2000f ) revealed traveling wave-like convection modes, 
which produce at the surface a combination of horizon- 
tal fiow velocities and pattern motions that have been 
associated with fiow properties observed in penumbrae. 
Recently substantial progress was made in "realistic" 
numerical simulations that include of the effects of par- 
tial ionization and radiative transfer. These models were 
first applied to sections of sunsp o ts (LS chiissler & Vdgfen 
2006; Hcincmann ct alj l2007t iRempel et al.l 12009^ 
Kitiashvili et al. 2009) and later to full sunspots 
([Rempel et ari l2009a: Rempelll20ldl ). 

Not all of the simplified models listed above con- 
tain a self-consistent description of the Evershed ef- 
fect, however, such a flow could be added as an addi- 
tional degree of freedom to most of them. The mod- 
els that include physical processes responsible for driv- 
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ing large-scale outflows are based on either stationary 
or dynamic flux tube models. In the case o f station- 
ary flux tube mod els (iMeve r fc Schmidt^ '1968*: 'Thomad 



1988; Deaenhardt 119891 11 991: Thomas & Montcsinos 
1993t iMontesinos fc Tho mas 1997), a pressure difference 
is imposed at the footpoints of the flux tube, which leads 
to siphon flows. In the dyn amic flux tube model of 
ISchlichenmaier et al.l (|1998al lbf) , fast outflows result from 
a combination of hot plasma rising at the inner footpoint 
and additional pressure driving that results from radia- 
tive cooling in the photosphere. 

In the more recent radiative MHD simulations, the 
penumbral fine structure is a byproduct of anisotropic 
overturning convection and the Evershed flow has been 
interpreted as the convective flow component in the 
direction of the magnetic field (Scharmer et al. 200^. 
A more detailed analysis of IRempel et al.l (|2P09a) con- 
cluded that convection in penumbra and quiet Sun dif- 
fer primarily in terms of anisotropy of the velocity field, 
while typical convective rms velocities and length scales 
of energy and mass transport are comparable. The con- 
vective nature of the penumbra is manifest in all ra- 
diative 3D simulations to date but it is still debated in 
the context of observational constraints. While there is 
evide nce for overturning convection in some investiga- 
tions (jZakharov et al.ll2008t lRimmelel[2008t iBharti et all 
[2010I) . many others primarily identify the upflow com- 
ponent of the Evershed flow in the inner and downflows 
in the outer penumbra with little e vidence fo r overturn- 
ing convection ( , Bellot Rubio et all [2005; .Ichimoto et al 



200faHlYa^z fc SchlichenmaieHl2009r iBellot Rubio et al 



2O10l) 



In this paper, we present a detailed a. nalysis of the 
model presented in IRempel et al.l (|2009a f) with special 
focus on the physical origin of large-scale outflows. After 
a brief description of the numerical model in Sect. [21 we 
focus first on the photospheric appearance of the penum- 
bra in Sect[3] In Sect. [3|-[5|we present a detailed analysis 
of the subsurface structure responsible for the driving of 
large-scale outflows. Sect. [6| analyzes the field geometry 
and connectivity in penumbral fiow channels and com- 
pares our results with findings from previous models, in 
particular models based on the flux tube picture. Sect. 
[7| analyzes deeper reaching outflow components beneath 
the penumbra that are not directly associated with the 
Evershed flow. The results are summarized and discussed 
in Sect. M 

2. NUMERICAL MODEL 

Our investigation is based on the simul ation of a pair 
of opposite polarity sunspots described in IRempel et al.l 
(l2009a). The simulation in a 98 x 49 Mm wide and 6.1 
Mm deep box was designed to study the formation and 
structure of penumbrae under a variety of different fleld 
strength and inclination angles. To this end the sim- 
ulation was initialized with a pair of opposite polarity 
sunspots each having a flux of 1.6 • 10^^ Mx each, but dif- 
ferent fleld strengths of about 3 and 4 kG, respectively. 
In order to focus on details of sunspot fine structure a 
rather high grid resolution of 3072 x 1536 x 384 (32 km 
horizontal and 16 km vertical) was used at the expense 
that this simulation could cover only a rat her short time 
span. While the original presentation in IRempel et al.l 
(|2009ari was based on a run of 1.5 h in high resolution 




X [Mm] 

Fig. 1. — Intensity image of the sunspot we analyze at about 
6 hours after initialization of the simulation run. We show here 
only one half of the horizontal extent of the computational do- 
main, which contains a pair of opposite polarity sunspots oriented 
in the x-direction. The panel is centered on the right sunspot in 
the simulation domain, which has a stronger magnetic field and is 
overall more coherent. Our analysis is focused on the penumbral 
region on the left side. The white lines indicate the region used for 
azimuthal averages. 

(3.5 h total) we have progressed the simulation in the 
meantime to 4 h in high resolution (6 h total). Dur- 
ing the extension of the simulation several aspects of the 
penumbral structure evolved. During the later stages 
of the simulation filaments became more radially aligned 
and the mean intensity profile shows in the inner penum- 
bra a constant value of about 0.7/© with a more steep 
drop toward the umbra. The overall properties of the 
Evershed flow did not show a signiflcant variation in the 
time frame covered by this simulation; we base our de- 
tailed analysis of the physical origin of the fiow pattern 
on the last hour of this simulation run (starting about 5 
hours after the initialization). 

We also emphasize for clarification that this simula- 
tion uses gray radiative transfer. When we refer in the 
following discussion to intensity, we mean the bolometric 
intensity, if we refer to optical depth, we mean the optical 
depth computed with the Rosseland mean opacity. 

3. PHOTOSPHERIC APPEARANCE 

3.1. Azimuthal averages in photosphere 

As described in detail in IRempel et al.l (j2009aD . the 
simulation domain contains a pair of opposite polarity 
sunspots, with the most extended penumbrae found in 
between the opposite polarity spots along the horizontal 
x-direction. The most coherent penumbra is found in the 
sunspot with the initially stronger field str ength of about 
4 kG ( see the spot on the right in Fig. 1 of lRempel et al.l 
l2009aD . We focus our detailed analysis on the latter, for 
which an intensity image is presented in Fig. [TJ We 
are here in particular interested in the extended penum- 
bra on the left side of the spot for which we highlighted 
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Fig. 2. — Azimuthal averages of quantities at the r = 1 level as function of the distance from the spot center. The dotted vertical lines 
indicate the radial positions also shown in Fig. [T] Quantities shown are a) intensity relative to quiet sun, b) radial velocity (solid) and 
vertical rms velocity (dashed), c) vertical field strength (solid) and radial field strength (dashed), and d) field inclination with respect to 
vertical. Here, the solid line indicates the average inclination angle, while the dashed line indicates the inclination computed from the 
average field. 



the sub-domain used for azimuthal averages in subse- 
quent figures. The dashed hues indicate i? = 10 Mm and 
i? = 20 Mm from the center of the spot. Fig. [2] summa- 
rizes properties at the r = 1 level, averaged azimuthally 
over the 90 degree wedge shown in Fig. [1] and about 1 
hour in time. The intensity normalized by the quiet Sun 
brightness 1q (panel a) shows a sharp increase from um- 
bra toward penumbra from about 0.15 to 0.7 /©. In the 
inner penumbra, the intensity stays constant on a plateau 
with about 0.7 Iq and increases then almost linear toward 
the edge of the penumbra where it reaches 0.95/© (due 
to the nearby opposite polarity spot in our simulation 
setup the intensity does not reach Iq). The plateau- like 
intensity profile formed during later stages of this simu- 
lat ion and was not presen t in the results reported earlier 
bv IRempel et all (|2009al) . The radial outflow (panel b) 
starts at about i? = 10 Mm, reaches its peak of about 
4 km near R = 15 Mm and drops off toward the outer 
edge of the penumbra, i? = 10 Mm corresponds to the 
position at which the average field inclination (displayed 
in panel d) angle exceeds 45 de grees, which was already 
found by IRempel et al.l ()2009aD as the critical value for 
the onset of large-scale outflows. The position of the peak 
velocity coincides with the position of maximum inclina- 
tion (about 70 degrees) in the middle of the penumbra. 
The inclination is defined here as arcsin(i3fl/|i3|). Due 
to the strong variation of inclination angle with azimuth, 
it makes a difference whether we compute the inclination 
locally and average in azimuth and time later or whether 
we base the computation on the averaged magnetic field 
presented in panel c). We show in panel d) both, the av- 



erage of the inclination (solid) and the inclination of the 
average field (dashed). The vertical rms velocity (panel 
b, dashed) increases steadily throughout the penumbra 
from a few 100 ms~^ at i? = 10 Mm to about 2kms~^ at 
the outer edge, which is the value corresponding to quiet 
Sun granulation. A value of about Ikms^^ is required 
near the inner edge of the penumbra (i? = 12 Mm) to 
maintain the penumbral brightness of 0.7 Iq. We find in 
this simulation an approximate relationship of the form 

3.2. Filamentation in photosphere 

Fig. [3] displays the filamentary fine structure seen at 
the T — 1 level in the penumbra. A filamentary structure 
is present in all quantities shown, however the strongest 
evidence is seen in intensity (panel a) , vertical magnetic 
field (panel c), inclination (panel d) and radial flow veloc- 
ity (panel e). Penumbral filaments show a strong reduc- 
tion of the vertical magnetic field strength, while the hor- 
izontal (radial) field component is moderately enhanced 
(panel b). The combination of the two leads to the 
strong variation of the inclination angle in the penum- 
bra. Strong radial outflows (panel e, red color indicates 
outflows) are seen in the almost horizontal flow channels, 
toward the outer end of flow channels the inclination an- 
gle exceeds 90 deg, indicating fleld returning back into 
the convection zone. Radial outflows with more than 
lOkms"'^ outflow velocity are indicated by solid con- 
tours. Most of the very fast outflows are found in the 
inner half of the penumbra, a few of them are associated 
with fast downflows in the outer penumbra. The verti- 
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Fig. 3. — Sunspot fine structure at the t = 1 level. Quantities shown are a) bolometric intensity, b) radial and c) vertical magnetic 
field, d) field inclination, e) radial and f) vertical flow velocities. A field inclination of Odeg corresponds to vertical field with the same 
polarity as the umbra, 90 deg to horizontal and 180 deg to vertical field with opposite polarity of the umbra. Radial outflows are displayed 
by red colors, solid contours indicate regions with more than lOkms"^ outflow velocity. Vertical upfiows are displayed by blue colors, solid 
contours indicate regions with more than Skms""*^ downflow velocity. 



cal velocity (panel f, blue colors indicate upfiows) shows 
strong up and downflows everywhere in the penumbra, 
strong radially aligned upfiows are preferentially found 
in the center of penumbral filaments. Solid contours in- 
dicate regions with more than 5kms~^ downflow veloc- 
ity. They are primarily found near the outer edge of the 
penumbra. We find downflow speeds of up to 15kms~^ 
near t — 1. Fast downf lows in opposite polarity region s 
have been observed by [Westendorp Plaza et al.l (120011 ): 
Idel Toro Iniesta et al.l (|2001l ). 

To clarify the relation between radial flow velocity, in- 
tensity and magnetic field strength in a statistical sense 
we present correlation coefRcients in Fig. Panel a) 
displays the correlation between intensity and radial ve- 
locity, panel c) the correlation between field strength and 
radial velocity. All correlations are computed based on 
the fluctuations of these quantities about their azimuthal 
mean. Intensity is correlated with outflows in the in- 
ner penumbra, but weakly anti-correlated further out- 
ward. The radial outflow is found in regions with re- 
duced field strength in the inner, but stronger field in 
the outer penumbra. S imilar correlations were found by 
(2007a) (see Fig. 3 therein) as well as 
(2005). In the panels on the right 



Ichimoto ct al. 
Schlichcnmaier et al. 

(b and d) we present additional correlations, which allow 



us to make a closer connection to the magnetoconvective 
structure of the penumbra. The radial dependence of the 
I — vji correlation is due to a decorrelation between ver- 
tical and radial velocity in the penumbra (panel b, blue) 
combined with a decorrelation of vertical velocity and 
intensity (panel b, red). While the latter remains posi- 
tive, a sign change is present in the former. The physical 
reason for the decorrelation between vertical and radial 
velocity is evident from the magnetoconvection pattern 
shown in Fig. [31 In the inner penumbra filaments are 
very narrow and the central upflow covers most of the fil- 
ament, leading to large positive correlation between the 
brighter upflow and radial outflow. In the outer penum- 
bra the patches of outflowing material become broader 
and several downflow lanes can be found within these 
patches, resulting in a reduction of the correlation. To- 
ward the outer edge of the penumbra stronger downflow 
patches are present, turning the correlation weakly neg- 
ative. Note that the I — Vz correlation stays low out- 
side i? = 20 Mm due to the proximity of an opposite 
polarity spot in our simulation setup. Fig. 2] panel d) 
presents additional correlations between radial velocity 
and vertical magnetic field (blue) as well as radial mag- 
netic field strength (red). While the former stays nega- 
tive throughout the penumbra, the correlation with the 
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Fig. 4. — a) Correlation between intensity and radial velocity ; b) correlation between vertical and radial velocity (blue) as well as vertical 
velocity and intensity (red); c) correlation between magnetic field strength and radial velocity; d) correlation between vertical field strength 
and radial velocity (blue) as well as radial field strength and radial velocity (red). 



radial field strength is positive. In the innermost penum- 
bra the strong reduction of \Bz\ dominates the picture 
and leads to an anti-correlation between vr and fur- 
ther out the contribution from the increased |J5/{| in the 
flow channels dominates and leads to a positive vr— \B\ 
correlation. The increase of the inclination angle found in 
the flow channels is a consequence of a strong reduction 
of Bz to almost zero, while Br is moderately enhanced 
in strength. This asymmetry is due to the fact that B^ 
benefits from a strong positive contribution of the induc- 
tion term {B ■ V)v due to the Evershed flow, while the 
corresponding term is negative for B^ due to the upward 
decreasing ve rtical velocity near r = 1 (see Sect. [5] for 
more detail). Ilchimoto et al.l (|2007aD found also a nega- 
tive VR — \B\ correlation in the inner penumbra with a 
trend of overall decreasing anti-correlation further out, 
however, a si gn change was n o t obse rved. The latter was 
proposed by iTritschler et al.l ()2007D and Ilchimoto et al.l 
|l008) based on observations of the net circular polariza- 
tion (NCP) in the outer penumbra at different viewing 
angles. 

The outflow velocity we find in the simulation is not 
stationary, we see flow variability that ranges from pe- 
riodic fluctuations on timescales of 5 — 10 minutes in 
the inner penumbra to quasi-periodic variations over a 
wider range of timescales starting from 10 — 20 min- 
utes in the center and outer penumbra. A flow vari- 
ability in the 15 — 40 minute range was a lso reported 
in the simulation of iKitiash vili et aT (2009) and associ- 
ated with_^v2ishe^^ 1994; Rimmele 
ll994HCabrera Solana et al.ll2007l ). It is conceivable that 
the periodic variations we find near the inner tip of fil- 
aments have a relation to twisting motions observed by 



Ilchimoto et al] ()2007b[ ) and lBharti eTall (|20Tol ). We fo- 
cus in this paper on the maintenance of the stationary 
flow component and base our analysis primarily on time 
and volume averages over sections of the penumbra. The 
non stationary flow component will be analyzed in a sep- 
arate publication. 

3.3. Mass and energy fluxes 

Fig. [5] panel a) displays filling factors of radial and ver- 
tical motions. While the upflow filling factor remains al- 
most constant around 0.4 to 0.5 from inner umbra toward 
the outer penumbra, the filling factor of outflows exceeds 
0.8 in the center of the penumbra. Panel b) shows the 
vertical rms velocity (black) , together with the mean ve- 
locity of upflow (blue) and downflow regions (red). The 
green line presents the mean vertical velocity averaged 
over regions with radial outflows (flow channels). The 
latter shows a weak average upflow of about 250 ms""'^ 
in the inner penumbra and a downflow reaching veloci- 
ties of more than 500 ms^^ toward the outer edge of the 
penumbra. 

The contributions from small- and large-scale flow 
components to mass and energy flux in the penumbra 
are presented in Fig. [51 In order to properly compare 
up- and downflow components we perform the analysis 
here on a constant height surface that is located about 
350 km beneath t = 1 in the quiet Sun (about half a 
Wilson depression downward). We decompose here the 
mass flux into positive and negative as well as azimuthal 
average components. Their contributions as function of 
radius are presented in panel a), where we show (mj") in 
blue and (mj) in red as well as (m^) in green (the latter 
is the sum of the former two). Here, (. . .) denotes the 
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Fig. 5. — a) Filling factors of radial outflows (red) and upfiows (blue), b) Vertical rms velocity (black), mean upflow (blue), mean 
downflow velocity (red), and moan vortical velocity in regions with outflows (green). 
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Fig. 6. — a) Mass flux in upflows (blue), downflows (rod), and azimuthal average of mass flux (green). In the inner penumbra up to 50% 
of the total upward directed mass flux is found in the large-scale flow component. The vertical mass flux is balanced within R < 20.9 Mm. 
In this region about 13% of the unsigned mass flux is found in the large-scale flow component, b) Decomposition of convective energy flux 
into normalized small-scale energy flux (blue) and large-scale energy flux (red). In green color is shown the normalized surface intensity 
for comparison. Integrated over the region with balanced vertical mass flux (R < 20.9 Mm), the contribution from the large-scale flow 
component is 12%, consistent with the mass flux contribution. 



azimuthal average and mf = {rriz ± \mz\)/2. While in 
the innermost penumbra up to about 50% of the mass 
flux is present in the azimuthal average component, this 
fraction drops steadily toward the center penumbra. In- 
tegrated over the region i? < 15 Mm about 1/3 of the 
total upward flowing mass is found in the azimuthal com- 
ponent, while the major fraction (2/3) is still overturn- 
ing laterally. The mass flux in the penumbra is balanced 
within R < 20.9 Mm. Integrated over this region the 
unsigned mass flux in the mean component constitutes 
about 13% of the total unsigned vertical mass flux in the 
penumbra. 

Evaluating the relative contributions from large-scale 
and small-scale convective motions to the total convec- 
tive energy flux in the penumbra requires an appropriate 
decomposition of the vertical mass flux. A separation 
just into azimuthal mean and the respective fluctuation 
would not be sufficient since the latter assumes that the 
large-scale flow is axisymmetric and equally considers fil- 
aments with higher temperature and the region in be- 
tween with lower temperature in the enthalpy flux. The 
consequence would be an underestimation of the over- 
all contribution from the large-scale flow (in the region 
B < 20.9 the net contribution would be —0.07X0). In- 
stead we construct the vertical mass flux of the laterally 
overturning flow component mf as follows: in regions 
with positive (mz), we reduce the amplitude of upflows 



such that they are in a mass flux balance with downflows, 
in regions with negative {rriz), we reduce the amplitude 
of downflows such that they are in balance with the up- 
ward directed mass flux. The large-scale mass flux is 



then given by 



TO, 



if. Unlike the decomposition 



into azimuthal mean and corresponding fluctuation, this 
procedure does not change the position of upflow and 
downflow regions for mf and compared to to^ . With 
H — (eint +p)/f? + f^/2 we can now compute the energy 

S / L S / L 

flux components Fz = {mz H) , which are displayed 
in Fig. [6] panel b). While matches the intensity 
profile very well in the outer penumbra, there is a clear 
deficit present in the inner penumbra. The deviations in 
the umbra are due to the fact that our horizontal slice is 
located above the r = 1 in the umbra and therefore the 
convective energy flux is zero. The large scale energy flux 
F^ has an amplitude of about +-/jq in the inner penum- 
bra and —Lq in the outer penumbra. The relevant quan- 
tity here is the net contribution after carefully balancing 
the upflow in the inner and downflow in the outer penum- 
bra. Integrating F^ over the region R < 20.9 Mm (in 
which the large-scale mass flux is balanced) leads to a net 
contribution of 12% to the total convective energy flux. 
The latter is very consistent with the relative mass flux 
contribution of 13% we found before. In an integral sense 
the large-scale flow contributes only a small fraction of 
the energy radiated away in a sunspot penumbra, but 
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Fig. 7. — Radial and depth dependence of radial flows in penum- 
bra. The radial flow velocity is normalized with the rms velocity 
of convective motions outside the sunspot (which is shown in Fig. 
[H] panel a). The dark solid lines indicate the r = 1 and t = 0.01 
levels. The strongest outflows reach velocities comparable to the 
rms velocity and are concentrated around r = 1. We see a tran- 
sition from outflow to inflow at about r = 0.01 in the inner and 
at T = 0.001 in the outer penumbra. In deeper layers, especially 
in the outer penumbra, outflows with speeds reaching 50% of the 
convective rms velocity are present. The vertical dotted lines indi- 
cate sub-regions we use for a more detailed analysis in the following 
discussion. 

locally the contribution can be larger. If we use the dif- 
ference between /Lq (blue) and ///© (green) in Fig. 
[6] as a rough estimate for the missing energy flux we iden- 
tify a contribution of up to 50% in the inner penumbra. 
Note that we avoided in the above discussion the asso- 
ciation between Evershcd flow and the large-scale flow 
component since there is no clear deflnition of what the 
former encompasses. If we associate the Evershed flow 
only with the horizontal flow pattern that corresponds to 
upflows in the inner and downflows in the outer penum- 
bra we would conclude that this flow pattern plays only 
a minor role in the penumbral energy transport. This 
deflnition would essentially correspond to the "sources" 
and "sin ks" of the Evershed flow t hat have been iden- 
tifled b y iRimmele fc Marinol ()2006D and llchimoto et al.l 
(|2007af) . Also note that the contribution from large-scale 
flows increases with depth as all intrinsic scales of con- 
vection increase with depth. Properly quantifying their 
contribution in deeper layers requires numerical simula- 
tions over longer timescales (since convective timescales 
increase with depth), which is beyond the scope of the 
current investigation. 

4. SUBSURFACE FLOW STRUCTURE AND UNDERLYING 
DRIVING FORCES 

4.1. Flow structure beneath penumbra 

Fig. [7] presents the subsurface outflow structure as 
function of depth and radial distance from the center of 
the spot. The depth is measured relative to the aver- 
age height of the t = 1 level in the quiet sun. We will 
use the same height scale in all of the following figures 
except Fig. [161 where we use the average t — 1 level 
in the penumbra as reference. Flow velocities are nor- 
malized by the rms velocity found outside the sunspot 
at the corresponding height level (see also Fig. [8]). We 
have chosen this normalization in order to compare flow 
fields found in the penumbra to convective flows found in 
almost undisturbed convection. We refer to this velocity 
reference in the following as v^^^. While the outflow ve- 



locity stays around 0.4 — 0.5 v'^^^ in the deeper layers, the 
near surface layers stand out with flow speeds exceeding 
uji'j^g. The two different scaling regimes of the outflow 
velocity found in the near surface layers (uppermost 500 
km) and the deeper part of the domain indicate already 
different physical driving mechanism at work, which we 
will analyze further in the following discussion. We ex- 
clude here the lower most 2 Mm of our domain which 
arc partially influenced by the bottom boundary condi- 
tion. The solid black lines indicate the average r = 1 and 
T = 0.01 levels. The radial outflow velocity peaks close 
to r = 1 and falls off rapidly with height. An outflow 
is present to about r = 0.01 in the inner and r = 10"'^ 
in the outer penumbra. The azimuthally averaged mass 
flux changes sign between r = 10"'^ and r = 10"''', 
since it puts more weight on the region above the more 
dense fllament channels with fast outflows. Overall the 
simulation indicates that radial outflows in the penum- 
bra are expected to be found in the deep photosphere, 
which is consistent wit h recent spectropolarimetric in- 
versions (IS chlichenmaicr eraL|[200l iBellot Rubio et al.l 
12006'; Borrero et al. 2008^, but not earlier work by 
Rimmclc (1995) and Stanchficld et al.l (|T997 ) where el- 
evated flow channels were inferred. Whether the inflow 
above r = 10"'^ could be related to the inverse Ever- 
shed flow observed in the Chromosphere (jPialetis et al.l 
19851) is currently a n open question, even though also 
Borrero et al.l (j2008f) found observational evidence for 
an inflow nea, r tempe rature minimum. As described in 
iRempel et all (|2009bD we switch for reasons of numeri- 
cal stability to an isothermal equation of state in regions 
with P = Pgas/Pmag < 10~^ and limit the Lorentz force 
such that the Alfven velocity does not exceed 60kms~"'^ 
to prevent stringent time step constraints. The latter 
two could possibly influence this flow pattern near the 
top boundary while the influence on flows in the photo- 
sphere is rather weak. Despite substantial velocities of 
a few kms~^ the associated mass and momentum flux 
is negligible compared to photospheric flows due to the 
sharp drop in density. 

The vertical dotted lines indicate 3 regions we refer to 
in the following analysis. We consider the region in be- 
tween i? = 10 and R — 12 Mm as inner penumbra, the 
region in between R = 12 and R = 18 Mm as center 
penumbra and i? = 18 and i? = 20 as outer penumbra. 
We have chosen the boundary for the inner penumbra 
based on the intensity proflle (Fig. [2|) that reaches a 
value typical for a penumbra of 0.7/0 at i? = 12 Mm. 
Since our outer penumbra might not be fully represen- 
tative for conditions in a "typical" outer penumbra due 
to the presence of a nearby opposite polarity spot with 
an Evershed flow in the opposite direction, we separated 
out regions with i? > 18 Mm. R = 18 Mm is also the dis- 
tance where most of the dominant filaments of the center 
penumbra end (see Fig. [T]and[3]). 

Fig. [5] compares rms velocities in the plage region sur- 
rounding the sunspot and the center penumbra. The top 
panels present the absolute rms velocities, the bottom 
panels relative to v'^^^. Blue indicates the vertical rms 
velocity, green and red the horizontal components (green 
is along the filaments in the case of the penumbra) . In the 
plage region (more or less undisturbed convection) about 
half of the kinetic energy is found in vertical motions. 
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Fig. 8. — Vertical profiles of rms velocities in the plage region surrounding the spot (a, c) and in the penumbra (b, d). The latter rms 
velocities are computed in the center penumbra in between R = 12 and R = 18 Mm. Panels a) and b) show the absolute values, while 
velocities in panels c) and d) are normalized by the rms velocity profile outside the sunspot near the edge of the computational domain. 
Black indicates the total, blue the vertical rms velocity, green and red the 2 horizontal components. In the case of the penumbra green is 
the component along the filaments. In the penumbra there are 2 distinct regimes. In the deep layers (< — 1 Mm) flows are anisotropic, but 
show a depth dependence similar to convective flows in the plage region (i.e., they scale proportional to In the near surface layers, 

this scaling is still present for flows turning over laterally (blue and red curve), while the flow component along filaments shows a much 
steeper increase with height. 



the other half equally distributed among the horizontal 
components. This scaling holds very well over the 3 or- 
ders of magnitude in pressure stratification shown here 
(panel c). In the penumbra the vertical and horizontal 
rms velocity perpendicular to the filaments show a simi- 
lar scaling and relative strength, but overall their ampli- 
tude is reduced to about 40% of the respective values in 
the plage region. The rms velocity along the filaments is 
strongly increased with respect to the vertical rms, indi- 
cating a strong degree of anisotropy. The rms velocity in 
the direction of filaments is proportional to in more 
than 1 Mm depth and shows a steep increase toward the 
photosphere. Overall the kinetic energy is reduced in the 
deeper layers, but doubled in the near surface layers com- 
pared to the plage region. The apparent excess of kinetic 
energy found in the Evershed fiow compared to the plage 
region is due to a vertical redistribution of kinetic energy 
combined with anisotropy of the flow. 

4.2. Underlying driving forces 

In order to investigate the physical processes that lead 
to the driving of large-scale outflows around sunspots 



we analyze the energy conversion terms in the kinetic 
energy equation. Starting from the momentum equation 
we derive the following energy balance (we drop the time 
derivatives since we are interested in time averages): 



V ■ {gg- Vp)^ + V ■ {jx B) 

pressure/buoyancy Lorcntz 

~gv- [{v ■ V)v\ + v_F,r^ = 

Acceleration Viscosity 



(1) 



Under the assumption of stationarity the acceleration 
term is identical to the negative divergence of the ki- 
netic energy flux, gUv'^ /2. A negative acceleration term 
implies positive divergence, i.e. the volume element is a 
source of kinetic energy. In Fig. |9] we compare the dif- 
ferent contributions to the energy equation for the plage 
region (panel a) and penumbra (panel b). On a qual- 
itative level there is a large degree of similarity: Pres- 
sure/buoyancy forces are the main driver, close to the 
surface most of that energy input is used up by accel- 
eration forces, the remainder is balanced in about equal 
parts by work against viscous and Lorentz forces. In 
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Fig. 9. — Comparison of energy conversion terms in kinetic energy equation between plage region (a) and center penumbra (b). The 
quantities shown are averaged in volume and about 1 hour in time. Black: work by pressure/buoyancy forces; red: work by Lorentz force; 
blue: work by acceleration forces; green: work by viscous forces. Overall the pressure/buoyancy forces are the primary driver of flows. 
Close to the surface a major fraction of the energy is deposited into acceleration work, viscous losses and work against Lorentz forces have 
about equal contributions. 
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Fig. 10. — Plage region: Energy conversion terms in kinetic energy equation. The quantities shown are averaged in volume and about 
1 hour in time. Black: work by pressure/buoyancy forces; red: work by Lorentz force; blue: work by acceleration forces; green: work by 
viscous forces. Solid lines indicate upflow, dashed lines downflow regions, a) Contribution from terms in the vertical and b) from terms 
in the horizontal direction. In the vertical momentum equation pressure/buoyancy terms are in balance with acceleration terms. Most 
pressure/buoyancy driving takes place in the downflow regions. In the horizontal momentum equations pressure driving is offset primarily 
by work against magnetic and viscous forces. 



the penumbra the total amount of energy input by pres- 
sure/buoyancy forces in the uppermost 1.5 Mm shown 
here is reduced to about 40% and more concentrated to- 
ward the photosphere. The reduction in energy input 
is consistent with the overah reduced kinetic energy in- 
tegrated over this depth range. It is also notable that 
the work against the Lorentz force is not substantially 
different from the plage region (relative to the respective 
pressure driving) despite the quite different field strength 
and field structure. 

To investigate the driving of fiows in plage and penum- 
bra further we split now terms into the contributions 
from different grid directions as well as fiow directions, 
i.e. we consider the following 18 terms: 

Pt^{vf[gg.~{^pU) (2) 
i± = («±(/xB),) (3) 
At=-{Qvf[{v-V)v],) (4) 



Here, i indicates either the Cartesian directions x, y, z 
in the case of the plage region or the cylindrical compo- 
nents i?, $, z in the case of the sunspot penumbra. With 
vf — (wi ± |wi|)/2 we denote negative and positive ve- 
locity components. Note that we compute all forces on 
the Cartesian grid and use the transformation to cylin- 
drical coordinates only to separate the directions along 
and perpendicular to filaments in our nearly axisymmet- 
ric penumbra fragment (i.e. we compute terms like VrF^ 

and instead of VxF^ and VyFy with F being any 

one of the forces). The explicit expression for the vis- 
cous force is rather complicated due to the non-linearity 
of the underlying artificial viscosity scheme. In the fol- 
lowing discussion we do not explicitly compute the vis- 
cous terms, but indicate their approximate magnitude 
by using the quantity = -(P^ + Lf + Af). We 
confirmed a close relationship for a few snapshots, for 
which we restarted the code and extracted all numerical 
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Fig. 11. — Center penumbra (averaged in between i? = 12 to ij = 18 Mm): Quantities shown are the same as in Fig. 1101 In the vertical 
momentum equation (panel a) pressure driving is balanced by by work against the Lorentz force. Contrary to the plage region the pressure 
driving takes place in upflow regions. Acceleration terms are negligible. In the radial direction (panel b) outflows are primarily driven by 
magnetic forces (in the uppermost few 100 km), which are in balance with acc eleration terms. In deeper layers pressure driving becomes 
more important, but remains weaker than in the plage region (see also Fig. I12II . 




Fig. 12. — Energy conversion terms for the radial direction in the center penumbra on different scale (panel a). In more than 300 km 
depth flows are primarily driven by pressure forces. Pressure driving favors outflows, the Lorentz force opposes inflows more than outflows. 
In less than 300 km depth, the Lorentz force is the dominant driver for outflows, while pressure driving favors inflows. Panel b) shows 
the balance for horizontal motions perpendicular to the fllaments (azimuthal direction). Work by pressure forces is in balance with work 
against magnetic forces. The dotted line indicates in both panels the role of pressure driving in the plage region for comparison. 



dissipation terms. 

Formally the energy conversion terms are power densi- 
ties (work per volume and time) . For the sake of making 
the text more readable we will refer to them in the fol- 
lowing discussion very often as " work done by / against . . . 
forces" instead of "work done by/against ... forces per 
volume and time" . Since the former is simply the latter 
multiplied by a unit volume element and time interval it 
has no further impact on the physical meaning of these 
terms. 

Fig. [TU] shows the energy conversion terms for ver- 
tical motions (panel a) and horizontal motions (panel 
b) in the plage region. Note that we only show one 
horizontal direction due to isotropy. In the vertical 
direction pressure/buoyancy driving is in balance with 
work done against acceleration forces. Most of the pres- 
sure/buoyancy driving takes place in downflows due to 
their overdense material that cannot be supported by the 
pressure gradient. Pressure/buoyancy driving in upflows 



is much weaker since they are very close to a hydrostatic 
balance. Close to z = the sign of pressure/buoyancy 
driving is changing in upflows as a consequence of the 
overshoot layer in the upper photosphere. Magnetic and 
viscous forces play only a minor role in the vertical di- 
rection. Horizontal flows are primarily driven by pres- 
sure forces. Most of the energy is absorbed by magnetic 
and viscous forces, only a small amount is balanced by 
horizontal acceleration in the uppermost 500 km of the 
convection zone. 

In comparison to the plage region, the center penum- 
bra shows distinct differences (Fig. [TT|) . Almost all pres- 
sure/buoyancy driving takes place in upflow regions: the 
presence of strong magnetic field causes in the near sur- 
face layers to a steepening of the pressure gradient, lead- 
ing to almost hydrostatic balance in downfiows and ex- 
cess pressure driving in the upfiows. This excess pressure 
driving in upflows is in balance with work against the 
Lorentz force while in contrast to the plage region verti- 
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Fig. 13. — Energy conversion terms as functions of radius and 
depth in the penumbra. We show the differences between the en- 
ergy conversion for outflow and inflow regions. Panel a) shows 



A„, panel b) Lj^ — L^, and panel c) — . The contour 



lines indicate regions with average outflows of more than 2 and 
4kms~^. 

cal acceleration of fluid does not play an important role 
at all. Most of the energy extracted by Lorentz forces in 
the vertical direction is deposited into outward accelera- 
tion of fluid along filaments. In that sense the Evershed 
flow is driven by vertical pressure forces in upfiows that 
are deflected into the horizontal direction through the 
Lorentz-force. This "deflection" is most efficient very 
close to the surface: integrated from 1 Mm (500, 250 
km) depth to the top boundary about 44% (64%, 97%) 
of the pressure driving in the vertical direction ends up 
as Lorentz force driving in the radial direction. This way 
kinetic energy that is normally deposited into the vertical 
direction (plage region) gets transferred directly into the 
horizontal direction and accounts for the high anisotropy 
seen in the penumbra. 

Also the role of horizontal pressure driving differs sub- 
stantially from the plage region, which can be seen in 
Fig. 112b ). where we show the energy conversion terms 



for the radial direction on a different scale. Pressure 
driving is dominant below 300 km depth, but the over- 
all magnitude remains smaller than in the plage region 
at comparable depth. While pressure driving shows a 
preference for outflows in more than 300 km depth, it 
prefers inflows further up. Pressure forces are the pri- 
mary cause for the deep flow component with velocities 
of about 0.5 we identifled in Fig. [3 but their overall 
role for the near surface flow is limited: integrated from 1 
Mm (500, 250 km) depth to the top boundary the contri- 
bution from P+ relative to L+ is 100% (38%, 9%). If we 
consider only the components of the driving that break 
the symmetry between in- and outflows, — and 
— i^, the corresponding values are 13% (1%, —9%). 
Note that most of the contribution to P^ — P^ comes 
from the region i? > 17 Mm, further inward P^ — P^ is 
close to zero (see also Fig. [T3| . 

However, pressure forces remain the dominant driver 
for flows perpendicular (^-direction) to filaments (Fig. 
[T^b). Here, pressure driving is offset by work against the 
Lorentz force, while both acceleration and viscous terms 
do not contribute substantially. 

The different role of pressure driving compared to the 
plage region is primarily a consequence of anisotropy in 
terms of radially elongated convection cells in the penum- 
bra: radial pressure gradients are reduced, while lateral 
pressure gradients are enhanced compared to isotropic 
granulation. In addition the steepening of the vertical 
pressure gradient (that leads to the shift of pressure driv- 
ing from down to upfiows) results in an overall reduction 
of pressure close to the photosphere in comparison to the 
ambient stratification. 

The clear association between Lorentz force driving 
and the near surface fiow pattern is evident from Fig. 1131 
Here, we present the quantities — A~^, — and 
P^ — Pfj as function of radius and depth. Strong nega- 
tive values of Ajj — indicate outward acceleration of 
fiuid. These regions are confined to a narrow layer near 
T ~ I. Here, the Lorentz force is the primary driver, 
pressure terms have weakly negative contributions (they 
favor infiows) . The peak of Lorentz force driving and ac- 
celeration is found in between i? = 12 and i? = 15 Mm. 

In deeper layers pressure terms are in approximate bal- 
ance with Lorentz force terms, resulting in only minor ac- 
celeration work despite their overall amplitude. Toward 
the outermost edge the Lorentz force is overcompensat- 
ing outward directed pressure driving resulting in decel- 
eration of radial outfiows (positive values of — A'^). 

For comparison with Fig. [TT]we present in Figs. [T4land 
[T5]the same quantities for the inner penumbra (from R — 
10 to i? = 12 Mm) and outer penumbra (from i? = 18 to 
i? = 20 Mm). 

In the inner penumbra we see a forcing pattern that is 
very similar to that we found for the center penumbra, in 
particular with respect to the near surface layers where 
the Evershed fiow is driven. Differences are present in 
deeper layers; here, the Lorentz force is also dominant 
in driving outflows and is actually driving these outflows 
against horizontal pressure forces. 

In the outer penumbra, we see a fundamentally dif- 
ferent situation (which is in part a consequence of the 
nearby opposite polarity spot and the resulting strongly 
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Fig. 14. — Inner penumbra (averaged in between ij = 10 to ij = 12 Mm): Quantities shown are the same as in Fig. 1101 
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Fig. 15. — Outer penumbra (averaged in between _R = 18 to _R = 20 Mm): Quantities shown are the same as in Fig. 1101 



magnetized downflow lane in between). Here, outward 
directed pressure forces dominate the picture entirely, 
however, they do not lead to a strong outward accelera- 
tion of fluid. To a large degree they are opposed by the 
horizontal Lorentz force and the energy is transferred to 
the vertical direction, where the Lorentz force becomes 
the major driver for downflows. The latter is due to the 
fact that magnetic field in the outer penumbra turns back 
downward. We see in the uppermost layers only a weak 
signal from the horizontal Lorentz force driving outflows 
- this is expected since we are in the region where the 
Evershed flow declines quickly to zero. 

We will further discuss the deeper reaching flow com- 
ponent in Sect. [T] 

5. MAGNETIC FILAMENT SUBSTRUCTURE 
RESPONSIBLE FOR DRIVING THE EVERSHED FLOW 

5.1. Simplified momentum balance 

After describing the driving forces behind the Ever- 
shed flow in detail in the previous section, we present 
here a simple model of the underlying thermal, magnetic 
and velocity structure and reduce the overall picture to 
the most relevant terms in the equations. In order to 
carve out the typical structure of the regions responsi- 



ble for driving the Evershed flow we select regions which 
have both upflows and outflows. In Fig. [TBk ) we present 
the mean magnetic field, flow, and thermal structure as 
function of depth obtained by averaging over all such re- 
gions horizontally between R = 12 and R — 18 Mm. 
While the vertical magnetic field (black, dashed) is con- 
stant at about 400 G, the radial component (black, solid) 
increases monotonically from about 900 G to 1.8 kG at 
the T = 1 level (vertical dotted line) and drops again in 
higher layers to about 1.6 kG. The steep increase just 
below r = 1 is essential for the Lorentz force component 
driving the outfiow as we will describe below. The ver- 
tical velocity (blue, dashed) increases monotonically to 
about 1 knis^^ just 100 km beneath r = 1, followed by a 
sharp decline to a few 100 ms~^ above r = 1. The radial 
flow velocity peaks right at t = 1, the maximum ampli- 
tude is about 5.5kms~^. The solid red line shows the 
mean temperature profile with the corresponding scale 
on the right. Panel b) presents volume averages of pres- 
sure driving in upflows (black) work against Lorentz force 
in upflows (red) and work by Lorentz force in outflows 
(blue) and work against acceleration forces (green) for 
the same region. The dotted black line shows contribu- 
tions from horizontal pressure gradients multiplied by a 
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Fig. 16. — a) Average vertical thermal, magnetic and flow proflles in regions responsible for driving the Evershed flowr. Magnetic fleld 
(black) and flow profiles (blue) are shown on the left scale, solid lines indicate the radial, dashed lines the vertical component. The 
temperature (red) is shown on the right scale, b) Pressure driving in upflows (black), work against Lorentz force in upfiows (red), work by 
Lorentz force in radial direction (blue) and work against acceleration forces in radial direction (green). The dashed lines show simplified 
expressions that are explained in the text. The dotted line indicates the contribution from horizontal pressure gradients multiplied by a 
factor of 10. Bottom panels: Contributions in induction equation from advection (black), field line stretching (red), fiow divergence (blue) 
and numerical diffusivity (green). Panel c) displays terms for the radial and panel d) for the vertical field component. Here, dashed and 
dotted lines refer again to simplified expressions explained in the text. Vertical dotted lines indicate the average r = 1 level. 

factor of 10. Averaged over the region shown (from -0.4 
to 0.2 Mm) they contribute about 10% to the total ac- 
celeration work. Work by the horizontal Lorentz force 
(blue) is 90% of the work by vertical pressure driving 
(black). The solid lines are based on all terms in the 
equations (see Eqs. ([U to dU ), the dashed lines are an 
approximation for Lorentz force and acceleration terms 
based only on volume averages of the following expres- 
sions: 



by a factor of 2, i.e. quantities such as QVRVzdzVR and 
QVRVzdzVR agree in general on a qualitative level for the 
region we selected to perform the averages. 

In the vertical direction we have essentially a magneto- 
hydrostatic balance involving the terms: 



d 



Svr 



-Q9 



(8) 



1 „ dBR 

-Vz — Br—— 

47r oz 



Lr = vr—Bz^— 
47r oz 



A, 



-QVrVz- 



dvR 
dz 



(5) 
(6) 

(7) 



This is evident from the opposing contributions of the 
terms —Vz {dzP + gg) and —VzBRdzBR/ (47r) in Fig. [TBb) 
(black and dashed red curve). The energy extracted 
by the Lorentz force in the vertical direction leads to 
a strong acceleration of an outflow in the radial direc- 
tion. Here, we have a balance between the Lorentz force 
and acceleration terms: 



The excellent agreement allows us to understand the 
driving mechanism behind the Evershed effect by con- 
sidering a reduced set of equations. Note that we could 
go even one step further with these simplifications, by 
expressing all terms through the mean quantities Br{z), 
Bz{z), vr{z), Vz{z), p{z), and g{z). Despite the fact 
that we are dealing with nonlinear terms of spatially 
highly inhomogeneous quantities, the agreement remains 
excellent except for the acceleration term that falls short 



At: ^ dz 



QVz 



dvR 
dz 



(9) 



The acceleration force results from the upward transport 
of plasma in a region with an upward increasing Evershed 
flow velocity. The work by vertical and radial Lorentz 
force components is in approximate balance, i.e. 
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leading to a simple relation between vertical and radial 
flow velocities of the form 

VzBr « vrBz ■ (11) 

The latter is the relation one would expect from a sim- 
ple " deflection" of vertical flows by an inclined magnetic 
field. 

Note that most of the acceleration of the outflow takes 
place about 100 km beneath the r = 1 level, while the 
outflow peaks right at r = 1. The latter is a consequence 
of the strong vertical upflow advecting accelerated fluid a 
few 100 km further upward. This upward advected fluid 
also overpowers the inward directed Lorentz force found 
right above r = 1 due to the sign change in OzBr. 

5.2. Induction equation 

Since the large positive value of dzBji right below t — 1 
plays an essential role in the acceleration process, we ana- 
lyze now how this magnetic field structure is maintained 
in the presence of strong vertical and radial flows. To 
this end we evaluate the different contributions in the 
induction equation: 

rifi 

— = -(Tj.y)B + {B ■V)v-B{V -v) . (12) 

Advcction Stretching Divergence 

The bottom panels of Fig. [T6]present these contributions 
(black: advection, red: stretching, blue: divergence) for 
the maintenance of the radial magnetic field structure in 
panel c) and vertical magnetic field structure in panel d) . 
As before solid lines show the full expressions, dashed- 
line approximations are described in the text below. For 
the radial field component (panel c) the dominant source 
is the stretching term in the induction equation. The 
major contribution to this term comes from the vertical 
shear profile of the Evershed flow, leading to an induction 
term BzdzVn (red dashed line). The remainder is due to 
horizontal stretching from terms like Budfivji. The peak 
of the stretching term (including all contributions) is lo- 
cated about 100 km beneath t — 1, where we also find 
the peak of the Lorentz force driving. This is not ex- 
actly where we find the peak of Bn, since there is an 
additional strong contribution from the vertical advec- 
tion that pushes strong radial field upward {—VzdzBn, 
black dashed line). The remainder is offset by the nega- 
tive contribution from the diverging convective motions 
(solid blue line). In the case of the vertical field (panel d), 
the role of the contributions from advection and stretch- 
ing are opposite. Here , advection is the primary mech- 
anism that maintains the field. The positive sign orig- 
inates primarily from horizontal advection terms (black 

dashed) with a dominant contribution from —vrOrBz 
due to the on average outward decreasing vertical field 
strength, but there are also positive contributions from 
vertical advection —VzdzBz- The dominant negative con- 
tribution to the stretching term is due to BzdzVz (red 
dashed line), which peaks close to r = 1 where dzVz is 
strongly negative. The remainder is offset again by the 
negative contribution from the diverging convective mo- 
tions (solid blue line). In both panels the green curve 
indicates the negative sum of these three terms, i.e. the 
amplitude of additional contributions from artificial nu- 



merical diffusivity. For both radial and vertical magnetic 
field the contributions from stretching, advection and di- 
vergence are in balance at the level of a few %. This indi- 
cates that the magnetic structure within the penumbral 
filaments in this simulations is not strongly affected on 
average by the unavoidable artificial magnetic diffusivity 
of the numerical scheme. Also the fact that our simula- 
tion contains almost field free umbral dots on scales even 
smaller than filaments sets strong constraints on the role 
artificial diffusivity plays. 



5.3. Filament cross section 

From Fig. [16] we deduce a vertical extent of about 
200 km for the region in which most of the energy con- 
version takes place. This value is obtained through an av- 
erage over all areas between R = 12 and R = 18 km that 
have upflows and outflows. Since the typical height vari- 
ation of the T = 1 level in the penumbra is about 200 km, 
this indicates that locally within individual filaments the 
fiow is driven in an even narrower boundary layer right 
beneath the r = 1 level. We illustrate this in Fig. [T71 
which shows magnetic field and velocity together with 
inclination and energy conversion by horizontal Lorentz 
force on a vertical cut through three developed and one 
just forming penumbral filament in the inner penumbra. 
It is evident that there is a narrow boundary layer form- 
ing along the t — 1 surface that is characterized by in- 
creased Br and reduced Bz, resulting in a strong increase 
of inclination. Lorentz force driving is concentrated to 
an equally thin layer just beneath t — 1, the resulting 
outfiow is broader and extends above t — 1. The latter 
is a consequence of the presence of overturning motions 
that transport and distribute accelerated fluid above and 
along the t = 1 surface. This redistribution does not re- 
quire additional acceleration work, since the associated 
kinetic energy flux is close to divergence free (the term 
for acceleration work is identical to the negative diver- 
gence of the kinetic energy flux under the assumption 
of stationarity) . In Fig. [T71 we highlighted a cross sec- 
tion in the inner penumbra, further out the fllamentary 
structure is less prominent. Nevertheless, we also see 
there a concentration of the energy conversion terms to 
very thin sheets beneath r = 1, while fast outflows are 
found mostly between r = 1 and r = 0.01. In a statis- 
tical (average) sense the differences between inner and 
center penumbra are small (compare Figs. 1111 and Il4 p . 
Also note that the concentration of energy conversion by 
Lorentz forces to thin sheets is typical for magnetocon- 
vection in a more general sense; however, the preferred 
location and quasi-steady maintenance of these regions 
near r = 1 is restricted to the penumbra. The mecha- 
nism responsible for the latter is explained in Sect. 15.41 



5.4. Formation of thin boundary layer 

In this section we illustrate the crucial role of the ver- 
tical advection terms by discussing a simplified model 
that captures the essential terms on a qualitative level 
to within a factor of two. We consider only those terms 
in the radial momentum and induction equations that 
have been identified as the dominant contributors in the 
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Fig. 17. — Vertical cross section of filaments in the inner penumbra. Displayed are a) radial and b) vertical field strength together with 
the resulting inclination in panel c). The bottom panels show d) radial, e) vertical velocity and f) the energy conversion by the horizontal 
Lorentz force along filaments. The two solid lines indicate the r = 1 and 0.01 levels. The enhancement of Bji combined with a reduction 
of Bz around r = 1 leads to a sharp increase of the inclination angle. The strong increase of the inclination angle is restricted to a very 
narrow boundary layer near t = 1, around which we also see the dominant contribution from the horizontal Lorentz force. The resulting 
outflow found around r = 1 is broader than the boundary layer, but also restricted to the deep photosphere. 
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Fig. 18. — Displayed are a) outflow and b) upflow velocities at t = 1 relative to |-Bn|/\/47rp and \Bz\/^/ATrg, respectively. Blue (white, 
red) colors indicate sub-Alfvenic (Alfvenic, super-Alfvenic) flow speeds. Throughout most of the penumbra the outflows stay close to 
0.8|Bfl|/v^. 
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previous discussion: 

dvB. 
dt 

dt 



dvR _ dBf 
^ dz 
OBr 



dz 



Ang dz 
dvR 
dz 



B 



(13) 
(14) 



Neglecting the advection terms, Eqs. (IT^ and dH]) lead 
to wave solutions of the form (assuming that B^ is nearly 
constant, which is at least true for the average shown in 
Fig. [HI): 



dt^ 



^d_ 

Aug dz-^ 



2 1 



ivR,BR)^0 



(15) 



In this case the profiles of vr and Br could not be main- 
tained in place and would spread out with the Alfven 
velocity corresponding to the vertical magnetic field com- 
ponen t, \Bz\/\/ATTg; see, for example. iVasil &: BrummeUl 
(|2009| ) for a discussion of the dynamics of magnetic 
shear layers. Since the Alfven velocity is of the order 
of 2kms~^ (using a mean value of B^ = 400 G and 
g = 3 ■ lO^^g/cm^ near the r = 1 level), the rather nar- 
row Evershed flow profile would broaden substantially 
within a few 100 sees of time. On the other hand, the 
inclusion of the advection terms allows for a stationary 
solution provided that Vz > \Bz\/\/AT:g. Since near r = 1 
the vertical field strength drops and upflows can reach 
locally up to 3kms~^, this condition can be met within 
the thin boundary layer in which most of the driving is 
taking place. With a sufficiently strong advection term 
the upflow counteracts the downward traveling Alfven 
wave, while the upward traveling wave is bound by the 
photosphere and transition to a low /3 regime. The quasi- 
steady maintenance of the shear layer despite the back- 
reaction of Lorentz forces seems to be at odds with Lenz's 
rule, however, we have to keep in mind that there is 
a steady flow of energy through the system ultimately 
driven by overturning convective motions. Indeed, the 
energy conversion by the vertical advection term in the 
induction equation, BrJ (Att) Vzd^BR, is identical to the 
energy extracted by the Lorentz force from convective 
motions in the vertical direction (Eq. [5]). 

To summarize, the most important feature responsible 
for driving the Evershed flow is a strong increase of Br 
just beneath the r = 1 level (combined with the presence 
of a vertical background field of a few 100 G). The steep 
gradient of Br is primarily maintained by the vertical 
shear profile of the Evershed fiow in combination with 
upward advection due to the strong upflow in the center 
of the filament. The vertical advection terms in the in- 
duction as well as momentum equation are essential for a 
quasi-stationary configuration with lifetimes far beyond 
the overturning timescale of convection. They ensure 
that the peak of vr and Br is maintained above the re- 
gion with the strongest Lorentz force driving as well as in- 
duction due to shear, which are proportional to dzBR and 
dzVR, respectively. The vertical confinement of the shear 
layer is guaranteed by the fact that locally the upfiow 
velocity can exceed \Bz\/ \/ATTg. Since the outflow veloc- 
ity is linked to the upflow velocity by the approximate 
relation vrBz ~ VzBr the expectation is that the result- 
ing outflow reaches a velocity of about vr = \Br\/ y^Airg. 
With Br = I... 2 kG and g = 3 ■ lO-'^gcm'^ the re- 



sulting velocities should he vr = 5 . . . lOkms^^, which 
is about the range we find for the velocity within flow 
channels. Fig. [18] displays outflow and upflow veloc- 
ities at T = 1 relative to \Br\/ y/Ang and \Bz\/ \/ATrg. 
Outflows are close to Alfvenic throughout the penum- 
bra (typically vr w 0.8 \BR\/y/A7rg), upflows are weakly 
super- Alfvenic. 

The fact that there is a clear threshold for the onset 
of this driving mechanism {vz > \Bz\/y/A'Kg near t — \) 
is a possible explanation for a more or less well defined 
inner edge of the penumbra, or related, a critical field 
inclination that needs to be exceeded (about 45 degrees 
in this simulation). Going further inward toward the 
umbra the vertical field becomes stronger while vertical 
velocities are reduced, which makes it harder to pass this 
threshold locally. Even if it would be passed the resulting 
radial flow velocities would be less due to the smaller 
value of Br. 

5.5. Observable consequences 

Unfortunately the " feature" responsible for driving the 
Evershed flow remains well hidden beneath the r = 1 
level. Even worse, if we compute the Lorentz force from 
the "visible" part of the magnetic field structure, the 
Lorentz force is inward directed due to the sign change of 
dzBR. The amplitude of the visible inward component 
above r = 1 is about a factor of 5 — 10 smaller than 
the strong outward directed component beneath t — 1, 
which is responsible for the outward acceleration. 

The observable parts of the magnetic and velocity field 
are an increase of vr and Br toward r = 1 and a strong 
vertical gradient of Vz around r = 1. The latter is only 
visible very deep in the photosphere. The combination 
of this three factors should lead to positive values of 
d|_B|/dr and d|t;ios|/dr for a variety of observation angles 
and therefore contribute to the net circular polarization 
observed in sunspot penumbrae. 

6. FIELD LINE STRUCTURE OF FILAMENTS 

Many simplified models of penum bral filaments and 
the or igi n of the E v ershed fl ow such as I Meyer fc Schmidtl 
(IT961: [Th OTmi (HOSl: IDegenhardtl (I1989L 11991: 



( 1997fl : iSchlichenmaier et al. ni998allbl ') are based on the 



iThomas fc Montesi iiosr(ll993 ); iMontesinos fc Thomai 



thin flux tube approximation. It is not clear from flrst 
principles whether penumbral filaments (fiow channels) 
can be identified with fiux tubes in a meaningful way. 
The latter assumes the existence of a well defined flux 
surface that encloses a flow channel and clearly sepa- 
rates fluid "inside" the channel from fluid "outside" and 
assumes further that there are well deflned "footpoints" 
when intersected with a horizontal plane somewhere be- 
neath the photosphere. 

The convective structure of the penumbra as presented 
in Sect. [3| puts already some limits on the usefulness 
of the flux tube concept, since overturning convective 
motions are mostly orthogonal to the flows assumed in 
the flux tube picture (except footpoints). These con- 
vective motions lead to a continuous mass, momentum 
and energy exchange while fluid is moving outward along 
penumbral flow channels. Integrated over the penum- 
bra this mass exchange is substantial, since we flnd only 
about 13% of the total unsigned vertical mass flux in the 
large scale flow component. 
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Fig. 19. — Field line connectivity and associated horizontal flow speeds in simulated penumbra. The color of the field lines indicates the 
radial flow velocity (the colors red, yellow, green, and blue correspond to velocities of < 0, 2, 4, and > 8kms~^, respectively). Filament 
1 indicates a peripheral umbral dot almost transitioning to a penumbral filament. Filaments 2-6 sample different radial position in the 
penumbra. The semi-transparent plane indicates a magnetogram near (average) r = 1. Smaller horizontal and vertical cross section indicate 
the regions from which we selected the seed points for the field line integration. 



In the following paragraphs we will discuss filaments 
on the basis of their field line structure and connectivity 
to allow for a better comparison with models that are 
based on the flux tube approximation. 

6.1. Field Line connectivity 

In Fig. [in] we present the magnetic field line connectiv- 
ity as well as radial outflow velocity. The filaments are 
representative for different radial positions in the simu- 
lated penumbra. The field line analysis presented here 
was performed u sing the VAPOR software package devel- 
oped at NCAR (IClvne fc RastJ [20051 : IClvne et al.l[200l 
(www.vapor.ucar.edu). The field lines are computed 
from a 15 minute average, which is about a character- 
istic wave crossing time along filaments. We have chosen 
the latter since in particular stationary flux tube models 
make only sense on timescales beyond that, but the fol- 
lowing conclusions are not affected by the averaging in a 
fundamental way. Filament 1 corresponds to a periph- 
eral umbral dot that almost transitions to a penumbral 
filament, filament 2 is representative for the inner, 3 and 
4 for the center and 5 and 6 for the outer penumbra. Seed 
points for filament 1 were chosen from a magnetogram 
at the umbral r = 1 level (regions with reduced field 
strength indicating peripheral umbral dot). Seed points 
for filaments 2-5 were chosen based on outflow velocities 
in the indicated vertical cross sections with more than 
5kms^^ (most of them are actually around 8kms^^). 
The seed points for filament 6 were chosen from the indi- 
cated horizontal plane based on regions with more than 
5kms~^ downflow speed. The color coding of the field 
lines indicates the radial outflow velocity (the colors red, 



yellow, green, and blue correspond to velocities of < 0, 
2, 4, and > 8kms~^, respectively). Going radially out- 
ward from filament 1 to filament 6 we see the following 
physical picture emerging. Near the umbra-penumbra 
boundary upflow plumes (similar to those forming um- 
bral dots in the center of the umbra) push mass upward 
along inclined field lines. The mass loading results in a 
small bend of the field line and the upflow is guided out- 
ward leading to outflow velocities of about a kms~^. Due 
to the strong almost vertical field the fiow is not power- 
ful enough to bend over field lines completely. This is 
consistent with a moderate enhancement of the field in- 
clination in peripheral umbral dots by some 10 — 20 deg 
that has been in ferred from spectropolari metric observa- 
tions (ISocas-Nav arro ct al. 2004; Riethmii ller et al.ll2008l: 
iSobotka fc J^ak 2009t lOrtiz et all l20]J) Going fur- 
ther outward (filament 2) field lines are bent over suf- 
ficiently to become horizontal for a distance of a few 
Mm. In the horizontal stretch we find outflows exceeding 
8kms~^ (blue color), nevertheless, the field lines remain 
connected to the top boundary. Near the outer edge of 
this filament, we see the formation of small dips right 
before the field lines connect back to the top boundary. 
The latter is a consequence of the mass flux decoupling on 
average from the field through the formation of U-loops 
and reconnection with deeper field lines that extend fur- 
ther out. In addition, downfiows present along the edge 
of filaments can transport field lines together with the 
mass flowing along them downward beneath the photo- 
sphere. Moving to the center penumbra (filaments 3 and 
4) the horizontal stretch with fast outflows is extended 
and field lines start to bend over and return beneath the 
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photosphere (filament 4). Here, niost of the mass un- 
loading from field lines happens through either the above 
mentioned U-loop formation or through field lines that 
bend over temporarily. Going further outward (filament 
5) the latter scenario happens most of the time resulting 
in a filament that follows the r = 1 level and returns 
back beneath the photosphere. Note that all of these 
filaments (2-5) have fast > 8 km outflows along their al- 
most horizontal stretch in the photosphere regardless of 
their connectivity further out. Filament 6 shows an ex- 
ample of field connectivity and fiows in proximity of one 
of the fast downflow patches in the outer penumbra (all 
of the selected field lines have more than 5 kms~ down- 
flow velocity near t = 1). Compared to the previous 
fllament most field lines host only outfiows in the 2 — 4 
km range with some faster fiows present in the ultimate 
proximity of the footpoint. Most field lines reach toward 
the higher photosphere and turn back beneath the pho- 
tosphere within 3 Mm. They do not show the extended 
horizontal stretches that host fast outflows in filaments 
2-5. 

6.2. Underlying physical picture 

Overall we do not see compelling evidence that the 
field line connectivity is linked to the presence of strong 
horizontal outflows. Filaments with more than 8kms~^ 
outflow speed can have any connectivity: in the in- 
ner penumbra fleld lines typically connect to the top 
boundary, in the outer penumbra field lines bend over 
and return beneath the photosphere. Looking at the 
smooth transition in the filament structure through- 
out the penumbra (filaments 2-5) strongly suggests that 
a similar process is responsible for driving outflows 
in all of them. Following up on the discussion in 
Sect, m and [5] outflows result from pressure driven up- 
flows that load field lines with mass and bend them 
over. A very similar s ituation was already described by 
Scharmer et al.l (1200 8^ based on the " slab" simulation of 
Heincman n et al.l (|2007) . The energetic signature of this 
process is a balance between pressure and Lorentz forces 
in upfiows and Lorentz and acceleration forces in the ra- 
dial direction. Horizontal pressure gradients do not enter 
the picture on average since the upflow cells are elongated 
in the radial direction, which strongly decreases the role 
of the radial pressure gradient. The elongation of fila- 
ments does not impact the Lorentz force since there only 
vertical field gradients matter (see Sect. [5]). The re- 
gion in which outfiows are driven is confined to a narrow 
boundary layer just beneath r = 1 as we discussed in 
Sect. [5] An interesting property of this driving mecha- 
nism is that no substantial acceleration work is present 
in horizontal stretches of the magnetic field, where we 
find most of the fast horizontal flows: the Lorentz force 
has no horizontal component there and horizontal pres- 
sure gradients do not contribute much. This is however 
no contradiction, since all of the fluid that appears at or 
above t — 1 has to pass through the narrow boundary 
layer with concentrated driving forces (see Fig. [T7| . The 
continuous vertical mass exchange along the fiow channel 
maintains an almost steady fiow despite the fact that no 
substantial driving forces exist above r = 1. 

The driving of outflows depends primarily on condi- 
tions in the upflow cell in the inner penumbra, the field 



line connectivity toward the outer penumbra is secondary 
and established as a consequence of the outflow (similarly 
also umbral dots form initially on field lines that might 
connect to a region several 100 Mm away, the field line 
connectivity changes as part of the process until over- 
turning convection is possible). The average field line 
connectivity found in the penumbra depends on the ra- 
dial position (combination of ambient field strength and 
inclination angle). Mass unloading happens primarily 
through U-loop formation and reconnection in the in- 
ner penumbra since fiows are not strong enough to com- 
pletely bend over field lines. In addition, entire field 
lines can be submerged by laterally overturning convec- 
tive motions. In the center and outer penumbra, contin- 
uous bending of the field lines increases the length of the 
almost horizontal part near t — 1 until mass can be un- 
loaded at the outer edge of the penumbra. This process 
happens periodically in the center and permanently over 
the life time of the filament in the outer penumbra. The 
fact that the flow speed in the flow channels does not 
show an increase when field lines bend over completely is 
a strong indication that the conditions of the outer foot- 
point are of secondary importance to the process. We 
cannot rule out additional contributions from unavoid- 
able numerical diffusivity allowing plasma to move across 
field lines; however, our analysis in Sect. [S] did not re- 
veal a very large contribution on average compared to 
the other terms in the induction equation. 

The above interpretation shares some sim ilarity with 
the "fallen fiux tube" concept of iWentzell ([1992) com- 
bined with the convect ive driving of outflows described 
here and previously by I Scharmer et all ()2008[ ). 

6.3. Implications for simplified models 

We address here only models that include driving pro- 
cesses for the Evershed flow, a more general discussion is 
presented in Sect. El 

The picture presented above shows substantial differ- 
ences to stationary siphon flow models that have been 
proposed to explain the Evershed flow. Those models 
assume that processes related to turbulent pumpin g near 
the o uter edge of the p enumbra (Montesinos fc Thomas! 
Il997t iBrummefl et al.l 120081) hold field lines down and 
establish the field line connectivity that allows then for 
siphon fiows due to pressure differences between the in- 
ner and outer footpoint. We see stronger evidence in 
our simulation for a fiow that is driven from within the 
penumbra regardless of the initial field line connectivity, 
even though siphon-like flow channels can result from 
this process in the outer penumbra (see, e.g., filament 5 
in Fig. [T9t . Since we find fast outfiows along horizontal 
stretches of field lines regardless of their field line con- 
nectivity we conjecture that siphon-like fiow channels in 
the outer penumbra are more a consequence of the fast 
outflow than its cause. 

This does not rule out additional contributions from 
processes as described bv iMontesinos fc Thomas! ()1997! ) 
in the outer penumbra. The fllament 6 we highlighted in 
Fig. [19] is a potential example for a siphon flow related to 
turbulent pumping near the edge of the penumbra. The 
footpoint in the outer penumbra is caused by a strong 
downdraft leading to an U-loop of field lines (see the up- 
ward returning fiux in the left corner of the indicated 
sub domain). The outer footpoint of filament 6 has as 
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Fig. 20. — Vertical (left) and radial (right) flow velocity at the r = 1 level (top) and 3.4 Mm depth (bottom). At the t = 1 level (3.4 Mm 
depth) the vertical flow velocity is saturated at 3kms~^ (0.8kms~^), the radial one at 6kms~^ (0.8kms~^). While vertical motions in 
the penumbra are subdued compared to the surrounding granulation, horizontal flows clearly stand out. Horizontal outflows are strongest 
in the x-direction, where the nearby opposite polarity spots impose a more horizontal fleld. In contrast to the flow fleld in the photosphere 
radial motions do not stand as much over the surrounding convection cells and are also much less dependent on the position (x vs. y 
direction). The presence of the sunspot leads to a ring-like arrangement of convection cells around the spot. In contrast to the more or less 
randomly arranged convection cells at larger distance from the spot this preferred arrangement leads to the generation of mean flows with 
an amplitude comparable to convective flows at the given height. 



a consequence fast downflows and low pressure (on aver- 
age more than lO^dynecm"^ lower than the regions most 
of the field lines connect to further inward). We see a 
bundle of arch like field lines extending a few 100 km 
above r = 1 and having outflow velocities mostly in the 
2 — 4kms~^ range (a few faster fiows are found in the 
proximity of the footpoint), which is consistent with the 
predictions of most stationary siphon models. However, 
the fiow velocities fall short of those present in filaments 
2-5, where fast outfiows are confined to almost horizontal 



stretches of the field lines in the deep photosphere over 
lengths of several Mm. When there are higher reaching 
arches present such as in filament 3 and 4, flow speeds are 
declining toward the highpoint, contrary to predictions 
from stationary siphon flow models. 

The concentration of driving forces to a very narrow 
boundary layer beneath r = 1 is not compatible with 
the acceleration of fluid along several Mm wide arches of 
field lines extending mostly above r = 1. 

Fast outflows in the deep photosphere are a 
natural outcome of the moving flux tube model 
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Fig. 21. — Vertical velocity and mass flux profile take at -R = 18 
Mm in the outer penumbra (cf. Fig. [T]!. The solid Une shows the 
radial flow velocity (scale on the left) and the dashed line the radial 
mass flux vjiq with the scale on the right. At a depth of 3 Mm 
the deep flow component transports about 10 times the mass of 
the Evershed flow due to the strong increase of mass density with 
depth. 

(jSchlichenmaier et al.lll998al lbl) in which processes sim- 
ilar to convective overshoot limit the vertical rise of 
plasma. The fast outflows found there have been at- 
tributed in part to a hot upflow near the inner footpoint 
that is magnetically deflected outward and in part to hor- 
izontal pressure gradients resulting from radiative cool- 
ing. The former has some similarity to the magnetic 
" deflection" we see in our magnetoconvection simulation 
leading to a process limited to a very narrow height range 
near r = 1, but we do not see significant contributions 
from horizontal pressure gradients in the radial direction. 
The primary difference is that in the simplified fiux tube 
picture there is only one inner footpoint present for the 
filament, while the driving process seen in our numerical 
simulation is approximately translation invariant along 
the filament, i.e. the entire filament is essentially a "foot- 
point" in which this process accelerates fluid. The trans- 
lation invariance also implies small contributions from 
horizontal pressure driving in the radial direction. 

7. DEEP FLOW COMPONENT 

In the previous section we focused on the driving of 
flows in the uppermost few 100 km of the penumbra. 
While these flows can be primarily explained through 
strongly enhanced Lorentz force driving, there is no 
strong contribution present in more than 400 km depth 
(cf. Fig. [T^ . The physical origin of this deeper reach- 
ing outflow becomes evident from Fig. [20l At the r = 1 
surface horizontal outflows stand out compared to typical 
granular flows at that height level, while vertical motions 
in the penumbra are subdued. Furthermore, strong out- 
flows are found preferentially along the x-direction where 
the nearby opposite polarity spots impose a more hori- 
zontal magnetic field. In about 3.4 Mm depth vertical 
and horizontal fiows do not stand out in terms of ampli- 
tude, but rather in terms of the overall fiow structure. In 
the periphery of the sunspot convection cells are arranged 
in a ring-like pattern in contrast to the random arrange- 
ment further away. This preferred arrangement leads to 
the appearance of large-scale mean flows (outflows away 
from the sunspot), which have amplitudes comparable to 



typical horizontal convective flows at the same depth, i.e. 
the flow amplitude should be a certain fraction of the rms 
velocity rather independent of depth as indicated in Fig. 
[71 While the appearance of an approximately axisym- 
metric mean flow is a consequence of geometric arrange- 
ment, the preferred outflow direction requires additional 
ingredients. A pure arrangement of convection cells in 
a ring-like fashion should lead to a pair of convection 
rolls, generating an inflow close to the spot and an out- 
flow further out. A preference for the outflow can be a 
consequence of the presence of strong magnetic field in 
the center, which inhibits motions converging toward the 
spot but has less influence on the diverging motions fur- 
ther out. In addition large-scale pressure gradients can 
lead to the preference of outflows. We see in this sim- 
ulation a combination of both. According to Fig. W2k k) 
radial flows in the center penumbra are driven primarily 
by pressure forces in more than 300 km depth. We see a 
preference for pressure driving of outflows, only a weak 
asymmetry is introduced by the Lorentz force, which op- 
poses inflows more strongly than outflows. Note that 
most of the pressure driving originates from the region 
in between i? = 16 and R—\% Mm, where the deep flow 
component gains speed. 

While the outflow close to r = 1 clearly dominates in 
terms of flow velocity, the deep reaching flow component 
transports significantly more mass than the shallow Ev- 
ershed flow. In Fig. [2ll we compare vertical profiles of 
(wfl) and {vrq) at the position i? = 18 Mm in the outer 
penumbra. At a depth of 3 Mm the mass flux is about 10 
times larger than the the mass flux of the Evershed flow 
in the photosphere. It is very likely that the deep flow 
component is related to large-scale moat flows observed 
around sunspots in the photosphere as well as deeper 
layers through helioseismology (see, e.g.. IGizon" fc BirchI 
pOO^ for a recent review and references therein). To 
clearly establish this relationship we need however a nu- 
merical simulation covering a substantially longer time 
span as well as depth range, which is beyond this inves- 
tigation. 

8. CONCLUSIONS 

We presented a detailed analysis of the recent numer- 
ical sunspot simulation by iRempel et al.l (|2009aD . Our 
investigation focused on properties of penumbral flne 
structure near the t = 1 level as well as the physical 
mechanisms behind the driving of large-scale outflows in 
sunspots at photospheric levels and beneath. Our main 
conclusions can be summarized as follows: 

• We find penumbral fine structure near t = 1, which 
is compatible with the observationally inferred pic- 
ture of an interlocking comb structure with fast 
> 8kms~^ Evershed flows along almost horizontal 
stretches of magnetic field. 

• Correlations between radial flow velocity, intensity, 
and field strength at the t — \ level show a good 
qualitative agreement with recent observations and 
are a direct consequence of convective energy trans- 
port in a sunspot penumbra. 

• The net contribution from large-scale fiows to the 
mass and convective energy fiux in the penumbra 
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is found to be about 12— 13%. Local contributions 
in the inner penumbra can reach 50%. Maintaining 
the penumbral brightness of 0.7Lq requires about 
1 kms~^ vertical rms velocity at the r = 1 level. 

• We find in the sunspot penumbra two flow compo- 
nents, which we separate according to their scaling 
with respect to the convective rms velocity out- 
side the penumbra While the deep compo- 
nent (more than 500 km beneath the photosphere) 
has an almost constant ratio to v'^^^ of about 50% 
rather independent of depth, the shallow compo- 
nent (uppermost 500 km) shows an increase toward 
the surface steeper than v^^^. While the latter is 
essentially the Evershed flow, the former is likely 
related to a deep reaching moat flow component. 

• The near surface flow component is almost exclu- 
sively driven through the horizontal component of 
the Lorentz force along filaments. The energy for 
maintaining this fiow is provided by vertical pres- 
sure forces in upflow regions, which are effectively 
deflected horizontally by an inclined magnetic field. 

• The Evershed fiow is driven in a thin boundary 
layer beneath t — 1. Essential ingredient is a 
strong vertical increase of Bji beneath t — 1 
combined with a moderate vertical average field 
strength of a few 100 G. The Evershed flow is 
strongly magnetized and reaches a peak velocity 
of about \B]i\/ y^Airg at t = 1. 

• Upward advection of momentum and magnetic 
field by overturning convective motions in the 
penumbra is crucial for maintaining the conditions 
under which a quasi-stationary Evershed fiow can 
be driven. 

• The deep reaching fiow component results from a 
preferred geometric alignment of convection cells 
in the periphery of the sunspot. Asymmetries in 
pressure and Lorentz forces lead to a dominance 
of the outflow component. The flow amplitude is 
about 50% of the convective rms velocity, rather 
independent of depth. 

• Flow channels cannot be easily identified with mag- 
netic flux tubes. The field line connectivity is 
changing between inner and outer penumbra and 
we see no compelling evidence that the field line 
connectivity plays a major role in determining the 
Evershed flow speed; on the contrary, fleld line con- 
nectivity is established primarily as a consequence 
of the outflow. 

A variety of simplifled as well as magnetoconvective 
models for the penumbra have been discussed to explain 
the Evershed effect. The majority of the simplified mod- 
els that have been used to describe the acceleration of 
horizontal fiows in the penumbra is based on stationary 
or dynamic flux tube models. 

We see strong limitations for the applicability of the 
thin flux tube approximation in the context of penum- 
bral flow channels. The continuous mass exchange along 
flow channels due to overturning convective motions is 



not part of the flux tube picture (it is essentially orthog- 
onal to the assumptions), but found to be substantial 
in the presented numerical simulation. In addition the 
changing field line connectivity along the flow channels 
from inner to outer penumbra does not allow for an easy 
identification with a flux tube, at least not for the whole 
length of it. It is in general difficult to find meaningful 
compact footpoints that are representative of the field in 
inner and outer penumbra at the same time. 

We see limitations for the applicability of sta- 
tionary flux tub e models such as Me yer fc Schmidtl 
(ilQfil: iThomasI (119881): IDegenhardtl' (119891. \'l99W 
iThomas fc MontesinosI (|1993[ ): iMontesinos fc Thomas! 
( 19971) as explanation for the flows in our simulation. 
Both, field line connectivity and the causality between 
outfiows and field line connectivity point toward pro- 
cesses in which the conditions in the outer footpoint (if 
it exists at all) are of minor influence on the flow pat- 
tern and outflows are mostly driven by convective mo- 
tions within the penumbra. In addition fast outflows are 
found preferentially in the deep photosphere along al- 
most horizontal stretches of field lines regardless of their 
connectivity, which is different from the situation de- 
scribed in most stationary siphon models to date. This 
does not rule out additional contributions in the outer 
penumbra from processes sim ilar to those described in 
IMontesinos fc Thomas! ()1997f ) , where turbulent pumping 
by convective motions at the periphery of the penumbra 
plays a crucial role in establishing the field geometry. 
Filament 6 in Fig. [19] is one possible example for such 
a configuration. Note that our simulation might not be 
fully representative of a "typical" outer penumbra due 
to the setup with a nearby opposite polarity spot, even 
though, the combination of strong horizontal field in be- 
tween both sunspots with a strong downfiow lane due 
to the converging Evershed flows from both sides has a 
tendency to enhance submergence of fleld by convective 
motions. 

The moving flux tube model of iSchlichenmaier et al.l 
(1998a,b) naturally produces outflows located in the deep 
photosphere and similarly to the situation in our simula- 
tion most driving is focused on the inner footpoint. The 
acceleration of outflows is attributed to a combination of 
" deflection" of hot upflows and horizontal pressure gradi- 
ents resulting from radiative cooling. In contrast to mov- 
ing flux tube models we see in our simulation strong lim- 
its on the overall role horizontal pressure driving plays for 
the acceleration of plasma in the radial direction, while 
the magnetic " deflection" of hot upflows is not restricted 
to the inner footpoint but found everywhere throughout 
the penumbra. 

Recently Th omas! (|2010l ) and also Priest (2010, private 
communication) have suggested that flows in the penum- 
bra could be described in terms of dynamical siphon 
models, in which pressure gradients are produced by 
MHD processes such as magnetoconvection and these 
drive time-dependent flows along the magnetic field, 
which then reacts by Lorentz forces to the presence of 
the fiow: in this view, field line connectivity plays only a 
secondary role and is more a consequence than a cause. 
While describing part of the picture, however, one also 
needs to understand the fluid motions that are responsi- 
ble for the filamentation and most of the energy trans- 
port. Furthermore, we showed that the magnetic driv- 
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ing of outflows originates from narrow boundary layers 
that form beneath t = 1 in regions where convective up- 
flows are present. These regions arise as a consequence 
of magnetoconvection and are not captured by flux tube 
models that do not include a filament sub-structure and 
o verturning conve ction. 

iGallowavl (|1975D presented a phenomenological model 
based on the " roll-convection" picture introduced by 
iDanielsonI (|1961[ ). which explains the Evershed flow as 
a consequence of unbalanced horizontal Lorentz force 
components: while the Lorentz force balances the gas 
pressure deficit of the sunspot on average, the filamen- 
tation of the penumbra implies a strong azimuthal vari- 
ation and therefore a violation of this balance locally 
that is responsible for driving the Evershed flow. It is 
argued further that the Evershed flow is located in dark 
filaments, where downward directed motions concentrate 
magnetic field and lead in return to an above average 
Lorentz force. While our magnetoconvection simulation 
certainly produces a Lorentz force that is strongly vary- 
ing in the azimu thal direction, our results disagree with 
IGallowavl ()1975D in detail as we find that most of the 
Evershed flow is driven in upflow regions. 

Several other penumbra models focus on the fine struc- 
ture of the penumbra, without necessarily addressing 
the driving of large-scale outflows. One of these mod- 
els which has gotten a lot of attention recently is th e 
"gappy" penumbra model by ISpruit &: Scharmen (j2006D : 
iScharmer fc SpruitI ()2006( ). In its original form this 
model proposed essentially field free upflow plumes em- 
bedded in an inclined background field. While this struc- 
ture certainly captures the essence of the convective pic- 
ture seen in the numerical simulations presented here, 
there has been a lot of discussion of how field free these 
"gaps" really are and to which extent the field strength 
seen in present simulations is affected by numerica,! diffu - 
sion (see, e.g., discussion in lNordlund fc Scharmeii[201Cl[ ). 
The physical explanation of the Evershed flow presented 
here implies the presence of strong Lorentz forces in the 
uppermost few 100 km beneath the t = 1 level in the 
penumbra, which requires the presence of strong 1 — 2 
kG magnetic field. In that sense our results are incom- 
patible with models that predict field free gaps at pho- 
tospheric levels! Furthermore we do not see evidence 
that the field structure is heavily infiuenced by numer- 
ical dissipation as the analysis presented in the bottom 
panels of Fig. |T6]reveals. Also it remains very controver- 
sial whether the wealth of spectropolarimetric observa- 
tions could be explained by an essentiall y unmagnetized 
Evershed flow (see, e.g.. iThomasI l2010l ). Nevertheless, 
the results presented here indicate an almost Salomo- 
nian solution to this discussion: Strong field is confined 
to a narrow boundary layer just beneath r = 1, while 
further down the field strength is substantially reduced 
compared to the ambient plasma (but still of the order 
of 1 kG). This scenario was also brought forward as a 
possible s o lution for this problem in a recent review by 
IScharmeri (|2009l ). 

Recently a variety of different magnetoconvection sim- 
ulatio ns with radiative transfe r such as | Heinema nn et al. 
(|2Ci07f) : iRempel et all ()2009b[ ): iKitiashvili etal.i (|2009il 
ha ve been used to model t he p enumbra. The models 
by iHeinemann et al.l (|2007D and IRempel et all ()2009b[ ) 
focused primarily on the transition from umbra toward 



inner penumbra, which corresponds roughly to the inner- 
most edge of the region we analyzed in this investigation. 
The energy conversion terms for that region are displayed 
in Fig. [T3I and do not show (except for the overall am- 
plitude) a fundamental difference to Fig. [TTJ From this 
we conclud e that the driving mecha nisms for horizontal 
outflow s in IHeinemann et al.l ()2007D and IRempel et al.l 
(|2009bD is essentially the same as discussed here. The 
overall picture we de sc ribed in Sect. [6] is similar to 



IScharmer et al.l (|2008l ). IKitiashvili et al 



2009D studied 



in an idealized setup the influence of field strength and 
inclination on large-scale flows and found a strong de- 
pendence of outflow speeds on the average inclination 
and fleld strength, which is consistent with the mecha- 
nism explained here. They also reported on temporal 
variations of the Evershed flow speed on timescales of 15 
to 40 minutes. 

Our analysis essen t ially reinforce s concl usions of 
IScharmer etahl (|2008D : IRempel et all (|2009aD that the 
penumbra is anisotropic magnetoconvection and that the 
Evershed flow can be understood as convective flow com- 
ponent in the direction of the magnetic fleld. The simi- 
larity between plage region and penumbra with respect 
to the different terms in the kinetic energy equation (Fig. 
[5]) is quite astonishing, a comparison between the depth 
proflles of rms veloci ties points tow a rd anis otropy as the 
main difference (see IRempel et al.l (|2009aD . supporting 
online material). Nevertheless, there are also notable 
differences which clearly differentiate the Evershed flow 
from horizontal flows in typical convection. While the 
latter is entirely pressure driven, the Evershed flow is al- 
most completely Lorentz force driven. Only flows that 
turn over laterally in penumbral filaments remain pres- 
sure driven. In addition pressure/buoyancy driving takes 
place primarily in upflow regions, in contrast to field free 
convection that is driven by top heavy downflow regions. 
An interesting new aspect pointed out in this paper is the 
confinement of the underlying driving mechanism to very 
narrow boundary layers that exist just beneath r = 1. 

For clarification we want point out that the driving of 
the Evershed fiow is achieved through the radial compo- 
nent of the Lorentz force while the total work done by 
Lorentz forces remains negative, i.e. the net effect is a 
sink for kinetic energy. The overall underlying energy 
source is convective instability, which enters the kinetic 
energy balance through pressure/buoyancy driving. The 
Lorentz force facilitates the energy exchange between the 
pressure driving in the vertical direction and the horizon- 
tal Evershed flow acceleration. A necessary condition for 
the latter is the shift of pressure/buoyancy driving from 
downflow to upflow regions in the penumbra we described 
above. 

In addition to the mechanism leading to the fast Ev- 
ershed flow in the upper most few 100 km of a sunspot 
penumbra we have also identifled a mechanism leading 
to the formation of a larger scale outflow in deeper lay- 
ers. In contrast to the Evershed flow the deeper flow 
scales proportional to the convective rms velocity (out- 
side the sunspot), the dominant radial outflow reaches 
typically amplitudes ~ 0.5i;rms- The main reason for 
this flow cell is a preferred circular alignment of convec- 
tion cells surrounding the sunspot. As a consequence the 
azimuthal average over this ring-like pattern of convec- 
tion cells does not vanish and leads to mean flow speeds 
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scaling proportional to t;rins- A preference for the out- 
flow results from a combination of pressure and Lorentz 
forces. While this flow does not stand out in terms of 
flow velocity as the Evershed flow, the radial mass flux is 
substantially larger, which make a connection with the 
large-scale moat flows observed around sunspots likely. 
In this simulation we do not see evidence for a converging 
collar flow that was found previously in 2D axisymmet- 
ricsirnulations (jHurlburt fc Rucklidgell2000tlBotha et al.l 
120061 120081 ). In a future publication we will investigate 
the subsurface structure of this flow component in deeper 
domains and evolution over timescales longer than those 
covered by the numerical simulation presented here. 

The deep flow component described in Sect. [7] should 
be in principle observable through local helioseismology. 
The clear prediction is here an outflow of plasma with an 
amplitude of about 50% of the convective rms velocity 
reaching downward several Mm beneath the penumbra. 
This result is in contrad iction with s ome recent h elioseis- 
mic inversions such as iZhao et al.l pOOll [2QToI ). which 
point toward an inflow in a depth range from 1.5 to 5 
Mm. One the other hand I Gizon et al.l ()2009( ) reported 
on an outflow over the uppermost 5 Mm. 

While most of the processes responsible for driving the 
Evershed flow are located beneath the r = 1 level, there 
are nevertheless several aspects of the magnetoconvective 
penumbra model presented here that can be constrained 
through observations. As presented in Sect. [3] most of 
the energy is transported in the penumbra by laterally 
overturning convective motions. We find a very tight re- 
lationship between intensity and vertical rms velocity of 
the form / cx a/ rms{T — !)■ From this follows that 
the vertical rms velocity at r = 1 in the penumbra with 
/ ~ 0.7/0 should be about half of the value found in 
the quiet sun, i.e. about 1 kms ~^. This value is consis- 
tent w ith the recent findings of iFranz fc Schlichenmaie"i) 
(|2009| ) , who computed from Hinode observations velocity 
distributions functions for both quiet Sun and penumbra 
(see Fig. 3 in their paper). The half width at half maxi- 
mum of the vertical velocity distribution function for the 
penumbra is about 500 ms""'^, while the same analysis 
results in Ikms""'^ for the quiet Sun, i.e. the latter falls 
short by about a factor of 2 compared to the value we 
find for the quiet Sun at r = 1 (due to a combination 
of limited observational resolution as well as the sharp 
decline of above t — 1). If we assume that the same 
shortfall applies also to t he penumbra, the vertica l veloc- 
ity structure reported in iFranz fc Schlichenmaied (|2009f ) 
is at least in a statistical sense fully consistent with the 
amount of overturning convection we see in the numer- 
ical simulation presented here. Other consequences of 
the magnetoconvective model are the sign changes in the 
I — vr and B — vn correlations presented in Fig. |4l 
They are consistent with the analysis of Hinode data pre- 
sented by Ichimoto et al. (2007_a) (see Fig. 3 therein) . A 
positive B — vr correlation in the outer penumbra was 
sugges ted by iTritschler et al.l ()2007D and llchimoto et al.l 
()2008() based on observations of the net circular polar- 
ization (NCP) at different viewing angles. As explained 
in Sect. [5] the observable consequences of the Evershed 
flow driving mechanism are a moderate increase of Br 
and a steep increase of vr toward t = 1 . Vz shows a very 



steep gradient in the deep photosphere. The peak veloc- 
ity of the Evershed flow in the deep photosphere should 
be around \BR\/y^4Trg. 

The fact that the simulated Evershed flow is a deep 
photospheric flow is a direct consequence of its convective 
origin. On the observational side the depth dependence 
of the E vershed fl o w is d ebat ed. While the inyestiga - 
tions by iRimmeld (|1995D and iStanchfield et al.l ()1997f ) 
point t oward flows in eleva t ed flow channels, recent 
work b ylSchlichenmaicr et al.' ('20041). iBellot Rubio et ahl 
(|2006( ). and Borrcro ct al. (2003) is in support of a 
flow in the deep photosphere declining with height. 
Another point heavily debated is the presence or 
absence of overturning convection in the penum- 
bra. Support for overturn ing conve ction is found by 
llchimoto et~an ( 2007b') : Zak harov et a l. (2008); Rimmcle 
('2008'): B harti et a l. (2010), whilelB cUot Rubio ct aL 
(2005); Ichimoto et al. (200731): IFranz fc Schlichenmaier 
(2003) ; iBellot Rubio et al. 11201 0[) see primarilv support 
for Evershed flow related upflows in the inner and down- 
flows in the outer penumbra - a flow pattern that ac- 
counts in our model only for a small fraction of the un- 
signed mass and energy flux integrated over the penum- 
bra. It appears that overcoming the discrepancy between 
the presence of overturning convection in MHD simula- 
tions and the lack of evidence in many high resolution 
observations is one of the biggest challenges both numer- 
ical models and observations will face in the future. It 
is unlikely that the absence of overturning convection is 
the solution to this discrepancy; a brightness of 0.7/© 
or more requires overturning mass flux at a level not 
much less than granulation. An other possible solution 
could be related to thin boundary layers, which are indi- 
cated but not well resolved in the simulation presented 
here. If quantities such as flow velocities and magnetic 
fleld change dramatically over short distances, moving a 
r-surface by a distance comparable to our grid spacing 
can make a dramatic differences for the visibility of such 
feature. This clearly indicates that the simulation pre- 
sented here can only be considered as a flrst step in that 
direction. 

A convergence study of the properties highlighted in 
this investigation covering the resolution range from 
96 X 32 to 16 X 12 km resolution (horizontal x verti- 
cal) is in progress. A preliminary analysis shows that 
most magnetoconvective properties of the penumbra are 
robust (qualitative agreement over the whole range in- 
vestigated, quantitative agreement from 48 x 24 km res- 
olution upward), while the photospheric appearance of 
sunspot flne structure improves substantially with reso- 
lution. The currently highest resolution case is presented 
in Rempel (2010). 
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